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Abstract 

Quark masses are of great prominence in high-energy physics. In this paper, we 
have studied the heavy meson systems via solving the Lippmann-Schwinger equation 
by using the Martin potential for heavy quark masses. We have also attempted to 
use Martin potential to find an acceptable mass spectrum for heavy quarkonia. We 
obtained this spectrum via minimal phenomenological model (M. Melles, Phys. Rev. 

D 62: 074019 2000). The mass spectra for bb and cc are calculated without taking 
into account the relativistic corrections and spin-dependent effects. The obtained 
mass spectra turn out to fit the experimental findings. By using the conventional 
spectrum, we extract the pole mass of heavy quarks and use it along with the 
relation between MS mass (modified minimal subtraction scheme) and the on-shell 
quark mass to find MS mass for these quarks. The obtained results for MS mass 
are in good agreement with corresponding values reported in the literature. 
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1 Introduction 


QCD has two important characteristics. First, due to asymptotic freedom, quark interac¬ 
tion in high energy is weak. Second, since quarks do not exit in isolation, the inter-quark 
force in low energies increases with quark distances [1]. Because of the high masses of 
such quarks as b , c, t in comparison to positronium, phenomenological potential models 
could be used for meson spectroscopy. Heavy quarkonium spectrum is a good means of 
investigating static interaction [2] and the mass of heavy quark is an essential requirement 
for exploiting this investigation. The quark masses are of prevalent interest among high 
energy physicists, however, it is impossible to measure directly the heavy quark masses 
because of the confinement of the quarks within the hadrons. The heavy quark masses 
could be defined indirectly through analyzing their effects on the hadrons [3]. 

The precious information about the quark-antiquark ( q-q ) interaction for ground and 
excited states of the mesons can be obtained by comparing the experimental data with the 
theoretical predictions [1]. In QCD Lagrangian, the quark mass parameters are not related 
to the quark masses which are used in the potential models. In these Lagrangians, there¬ 
fore, a more common definition of mass which depends on the renormalization scheme is 
used. Two different definitions of mass are often used in the renormalization of the QCD 
Lagrangian; MS mass and pole mass. The pole mass and the MS mass are renormalized 
quark masses in the on-shell renormalization scheme and in the modied minimal sub¬ 
traction scheme, respectively. The mass definition that we use depends on the physical 
situation [4], MS scheme is the most popular renormalization scheme for QCD perturba¬ 
tion. Non-perturbative effects at high energies or short distances are of little significance. 
The most QCD calculations are performed by using MS scheme which is applicable for 
both light and heavy quarks. This scheme is also of use in estimating the quark masses 
in Standard Model. For all the above, MS scheme is considered as a perfect scheme for 
studying the physical properties of quarks. 


The paper is organized as follow. In Sect. 2 we explain MS and pole mass and the 
relation between these two masses. In Sect. 3 we present the non-relativistic Lippmann- 
Schwinger integral equation for two quark bound state by using a local interaction. In 
Sect. 4 we have introduced a potential model for heavy quark systems and we have 
presented our numerical results for the mass spectrum and A4S mass for these systems. 
Finally a conclusion is given in the Sect. 5. 
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2 MS and Pole Mass 


In this paper we study the heavy quarks. For a heavy quark, m 3> A qcd [5], where A qcd 
is the scale of strong interaction in QCD, and its magnitude is approximately 200 MeV. 
We list the bare heavy quark masses used in this work in the Table 1. 


Table 1: The bare heavy quark masses. 


quark 

mass (GeV) 

c 

1.51 

b 

4.88 

t 

174 


One of the results of analyzing the heavy quark spectrum is that the position of the 
pole in the propagator is considered as the pole mass [4], The difference between the 
quark mass and the pole mass is considered as an imaginary part (which is a multiple of 
A qcd) [6]: 


m q — rn P oi Am, 

5m = iSm — i Im(m po i e ), 
Im(m po i e ) = const, x A QC d- 


The relation between these two masses is: 


rripoie = nriM S [ 1 + M(a s )\, 


where: 


M(a s ) = £ P n 


-n+1 


( 1 ) 

( 2 ) 

(3) 


(4) 


(5) 


ii=0 


where P n is a function of rq (rq is the quark flavor) and a s is the strong interaction 
constant. For three-loop order, the relation is[2, 7, 8 , 9]: 


m Q = m Q (m Q ) 


1 + 4 a s {m Q ) + / q^m . 

3 7T V 7r 


q) 


+ 6 


a Am. 


71 


q) 


( 6 ) 


where tuq is the pole mass, and uiQ^niQ) is the running mass in MS scheme and £2 
and £3 are a function of rq which are given explicitly in Ref. [10] as: 

£2 = 13.44 — 1.041 rq, 

£3 = 194(5) — 27.0(7) rq + 0.65 nf. 
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(7) 

( 8 ) 



If m po i e is substituted in terms of MS mass, then the static energy of the quark- 
antiquark system could be [11]: 

Etot ( r ) = 2 mpoie + Vqcd( r)- (9) 

Of course, the inter-quark force could be used instead of the total energy because the 
perturbative extension of F(r) is far more convergent than the potential [12]: 

F M = ~Yr EUr) = ~Tr VQCD{r) ‘ (10) 


3 Lippmann-Schwinger Equation for Two-Body Bound 
State 


The bound state of two particles which interact by potential V is described by homoge¬ 
neous Lippmann-Schwinger equation: 


= G,V\ 


( 11 ) 


G° = E * Hq is the free two-body propagator, where E b is the binding energy of two-quark 
bound state and H 0 is free Hamiltonian. In configuration space the equation (11) turns 
out as: 

Pjf r e -*s/m\E b \ |r—r'| 

0(r) = -m^-j d 3 r' -j- y - V(r') / ip( r'), (12) 

where m is the average mass of bare quarks. The compact form of the equation (12) is: 

r oo 

ip(r) — / dr'M(r,r')i() (r'), (13) 

Jo 

where: _ _ 

/ I — y/rn\Eb\ Vr 2 +r' 2 — 2rr'x' 

dx> — / 2 a o = r = T r ' 2 V ( r 0- (I 4 ) 

-1 yr + r — Zrr x 

The eigenvalue integral equation (13) can be written schematically as: 


K(E) \i>) = A (E) 


(15) 


where K(E ) stands for the kernel of eigenvalue equation, which is energy-dependent. A (E) 
is the eigenvalue, and 1-0) denotes the two-body wave function. For a physical solution of 
the eigenvalue equation 15 and to obtain the binding energy of the system, i.e. E = E b , 
one should solve this equation to obtain the eigenvalue A = 1. The integral equation 13 
can be solved by iteration or direct methods. Since there is not any shifted argument in 
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the wave function amplitude, i.e. ip(r'), it would be easier to solve this integral equation 
directly, without using any iteration procedure. In order to solve the integral equation 
13, we should first discretize the continuous configuration and angle variables [13]. To 
this aim, for discretization of both configuration and angle variables we use the linear 
mappings by using the Gaussian quadrature grid points, where r u r' and x' k are the mesh 
points corresponding to r, r' and x' variables with number of mesh points N r , N r i and 
N x i, correspondingly. By these considerations, the eigenvalue equation 13 can be written 
as: 

N r r 

(16) 

3 =1 

where the matrix elements of the kernel M (ry r'-) can be obtained from equation 14 as: 


N x' „ - y/rn | E b | y/rf +r' 2 -:2r; r'. x' k 

M(rvr') = -v^mvr 3 / 2 £ W x' k - , 

k= i \jrl + rf - 2r i r' j x' k 


r'?V(r'), 


(17) 


where W r > and W x > k are the point weights of Gauss-Legendre polynomials. In the last step 
by diagonalization of the kernel of eigenvalue equation: 


K{r i y j ) = W r . j M{r i ,r' j ), (18) 

one can obtain the binding energy of two-quark system for each energy level. 


4 Potential and Mass Spectrum 

Many potential models have been proposed for inter-quark interaction. Some of them 
suit for light hadrons, and some other are suitable for heavy hadrons. For example, the 
models proposed by A. De Rujula et al. are not used to describe the light-quark systems 
[14], Martin models apply well to nonrelativistic calculations of heavy mesons [15] and 
Cornel potential is appropriate for heavy and light hadron systems [16, 17]. 

In this paper, we have used the Martin potential (low-power potential) to study the 
mass spectrum of heavy mesons such as bb and cc systems [15]. The functional form of 
the Martin potential is given by: 

V (r) = b m + a m (c m r) ai , (19) 

where the parameters of potential are: 

b m — —8.064 GeV, a m — 6.898 GeV, c m = lGeV. (20) 
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Accordingly, in the first step, i.e. estimating the binding energy for bottomonium and 


charmonium, has turned out to submit reasonable findings based on the mass formula. We 
have shown our results for the mass spectrum of bottomonium and charmonium systems 
in Tables 2 and 3 in comparison to other theoretical results and experimental data. 

Table 2: Calculated bottomonium (bb) mass spectrum (in unit GeV) for Martin potential 
and for spin=l. 

state 

Ref. [18] 

our work 

experiment [19] 

IS 

9.46 

9.45 

9.460 

2S 

10.03 

10.051 

10.023 

3S 

10.36 

10.32 

10.355 

IP 

9.90 

9.89 

9.9 

2P 

10.26 

10.249 

10.26 

ID 

10.15 

10.169 

10.161 

Table 3: Calculated charmonium ( 

[cc) mass spectrum ( 

in unit GeV) for Martin potential 

and for spin=l. 




state 

Ref. [18] 

our work 

experiment [19] 

IS 

3.097 

3.07 

3.068 

2S 

3.69 

3.75 

3.672 

3S 

4.78 

4.09 

4.040 

IP 

3.528 

3.528 

3.525 

2P 

3.944 

3.963 

— 

ID 

3.806 

3.803 

3.779 


In order to be able to calculate the running mass in MS scheme for heavy quarks from 
equation 6, we need to extract the pole mass m po i e from the calculated mass spectra, our 
results for the pole masses are as: 


m b { 2 - loop) ~ 4.72 GeV, 
m c (2 — loop ) ~ 1.4 GeV. 

By using the following values for coefficients a s and £ 2) which are taken from Ref. [2], 
for b for c quarks: 

a s {fh b ) ~ 0.217 

6 ^ -9.27 (22) 
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(23) 


a s (m c ) ~ 0.31 
6 - -10.31 

and by considering the mass poles given in equation 21, we can obtain the MS masses, 
which are given in Table 4. The comparison of onr numerical results with the values 
reported by the Particle Data Group (PDG) [20, 21], indicates that the two sets of results 
are in good agreement . 


Table 4: The running mass in MS scheme for heavy quarks. 


meson 

Refs. [20, 21] 

our work 

rrTb(rnb)( 2 — loop ) 

4.19 

~ 4.33 

m7(r?r7)( 2 — loop ) 

1.23 

~ 1.35 

rn7(rn7)(2 — loop) 

166 

~ 166.36 


The top quark mass is an important parameter in standard model at high energies [22]. 
Top quark has a width of T = 1.5 cm, which is larger than A ocd- Because of this short 
lifetime, top quark decays very quickly. MS mass, which is a short-distance mass, could 
be determined very accurately [21], For t quark, we follow the same procedure. However, 
a mass spectrum does not exist for this unstable meson by having a short lifetime [21]. 

We have solved the Lippmann-Schwinger equation for this meson by using the Martin 
potential and we have shown our numerical results for tt binding energy for both spin 0 
and 1 states in table (5). 


Table 5: tt binding energy for both s=0 and 1 states. The parameter A denotes the spin 
coupling coefficient. The stability interval [13] is related to the configuration cutoff which 
has been considered in numerical solution of the eigenvalue integral equation 12. 


spin 

A[AdeV 3 /c 6 } x 10 10 //r 2 

BE [MeV] 

stability interval [fm 1 

S=1 

-449.01 

-11.296 

[0.3, 0.4], [0.44,1.2] 

s=0 

45.759 

-36.948 

[0.3, 0.4], [0.43,1.1] 


By using the following values for the coefficients a s and £2 [23]: 

a s (m t ) = 0.1085 

£2 - -10.31 (24) 
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we obtain the MS mass for t quark which is given in Table 4 and is in good agreement 
with corresponding values reported in Refs. [20, 21]. 

We would like to add the comment that the integral equation 12 is singular for the 
confining potentials, and consequently the calculated energy eigenvalues would not be in 
agreement with the exact analytic binding energies. To overcome this problem one can 
use a regularized form of the confining potentials to remove the singularity of the kernel. 
To this aim one can keep the divergent part of the potential constant after exceeding a 
certain distance, which creates an artificial barrier. So, we have fixed the potential in the 
stability interval and our numerical calculations show that the physical eigenvalue A = 1 
is quite stable in this region for configuration cutoffs. 

5 Conclusion 

We have used the Martin potential to study the heavy quark systems by solving the ho¬ 
mogeneous Lippmann-Schwinger equation without considering the relativistic corrections 
and also the spin-dependent effects. We have obtained the mass spectrum for the heavy 
mesons and our numerical results are in good agreement with other theoretical results 
and experimental data as well. We have calculated the pole mass by fitting the mass 
spectrum of b and c quarks to other references. We have specified the running mass MS 
for these systems. Our obtained masses are close to the findings reported by the Particle 
Data Group. For t quark, we have verified the stability interval of regularized form of the 
Martin potential. 
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